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^ ■ Abstract 

The adsorption of highly oppositely charged flexible polyelectrolytes on a charged sphere is in- 
vestigated by means of Monte Carlo simulations in a fashion which resembles the layer-by-layer 
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Q ' deposition technique introduced by Decher. Electroneutrality is insured at each step by the pres- 
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ence of monovalent counterions (anions and cations). We study in detail the structure of the 
equilibrium complex. Our investigations of the first few layer formations strongly suggest that 
multilayering on a charged colloidal sphere is not possible as an equilibrium process with purely 
electrostatic interactions. We especially focus on the influence of specific (non-electrostatic) short 
range attractive interactions (e.g., Van der Waals) on the stability of the multilayers. 
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I. INTRODUCTION 



Polyelectrolyte multilayer thin films are made of alternating layers of polycations (PCs) 
and polyanions (PAs). The so-called layer by layer method, first introduced in planar geom- 
etry by Decher, consists in a successive adsorption of the polyions onto a charged surface and 
has proved to be extremely efficient 2|. Due to the many potential technological applica- 
tions such as biosensing [3], catalysis optical devices |5] etc., this process is nowadays 
widely used. Various techniques are employed to control the po 



such as optical [6 
NMR techniques 



m and neutron 



reflectometry, AFM 
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10|, UV spectroscopy 
and others. Some experiments (see e.g., Ref. were devoted to 



the basic mechanisms governing polyelectrolyte multilayering on planar mica-surfaces where 
especially, the role of surface charge overcompensation was pointed out. 

Another very interesting application is provided by the polyelectrolyte coating of spher- 



ical metallic nanoparticles 



15l |. This process can modify in a well controlled way the 



physico-chemical surface properties of the colloidal particle. Despite of the huge amount of 
experimental works, the detailed understanding of the multilayering process is still rather 
unclear, especially for a charged colloidal sphere. Hence the study of polyelectrolyte multi- 
layering is motivated by both experimental and theoretical interests. 

On the theoretical side, the literature on this subject is rather poor. Based on 
Debye-Hiickel approximations for the electrostatic interactions and including lateral cor- 
relations by considering different typical semiflexible polyelectrolyte-layer structures, Netz 
and Joannvjiol found a remarkable stability of the periodic structure of the multilayers in 
planar geometry. For weakly charged flexible polyelectrolytes at high ionic strength qualita- 
;ive agreements between theory 17|, based on scaling laws, and experimental observations 
^, Q| have been provided. The driving force of all these multilayering processes is of electro- 
static origin and it is based on an overcharging mechanism, where the flrst layer overcharges 
the macroion and the subsequent layers overcharge the layers underneath. However, the 
role of non-electrostatic interactions though pointed out in Ref. 0, Q] is not clear. In 
particular, it is still open whether the layer build up is an equilibrium- or out of equilibrium 
process, which resembles more a succession of dynamically trapped states. Therefore we 
do not know whether or not the complex polyelectrolyte is in thermodynamical equilibrium. 
This point has also been emphasized in a recent experimental work on planar multilayers 
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0] where considerable kinetic effects were reported. So far, there are nor analytical results 
neither simulation data for multilayering formation onto charged spheres. 

The goal of this paper is to study the underlying physics involved in the polyelectrolyte 
multilayering onto a charged colloidal sphere by means of MC simulations. The paper is 
organized as follows: Sec. IHlis devoted to the description of our MC simulation method. The 
relevant target quantities are specified in Sec. IIIIl The single chain adsorption is studied 
in Sec. IIVI and the polyelectrolyte bilayering in Sec. Then the multilayering process is 
investigated in Sec. IVII The case of short polyelectrolyte chains is considered in Sec. I VIII 
Finally, Sec. IVII II contains some brief concluding remarks. 



II. SIMULATION METHOD 



The setup of the system under consideration is very similar to those recently investigated 
by means of molecular dynamics simulations [2^ 0| . Within the framework of the primitive 
model we consider one charged colloidal sphere characterized by a radius a (= 4.5o") and a 
bare charge Qm = —Zm^ (where e is the elementary charge and Zm = 40) surrounded by 
Zm neutralizing monovalent {Zc = 1) counterions and an implicit solvent (water) of relative 
dielectric permittivity ~ 80. In the remaining of the paper, we will refer to the term 
macroion as the charged colloidal sphere. Polyelectrolyte chains (A^+ PCs and PAs) are 
made up of monovalent monomers [Zm = 1) of diameter a. For the sake of simplicity, 
we only consider here symmetrical complexes where PC and PA chains have the same length 
and carry the same charge in absolute value. To each charged PC or PA we also add 
small monovalent {Zc = 1) counterions (anions and cations countering the charge of the 
polyelectrolytes) of diameter a, hence always a charge neutral entity was added. Thereby 
all the microions have the same valence Z = Zc = Zm = 1 as well as the same diameter 
cr. Added salt of course would even weaken the effects observed and would be especially 
important for the case of an adsorption interaction between macroion and polyelectrolyte. 

All these particles making up the system are confined in an impermeable spherical cell of 
radius R = 60a. The spherical macroion is held fixed and located at the center of the cell. 
To avoid the appearance of image charges ||22|], we assume that the macroion has the same 
dielectric constant as the solvent. 

Standard canonical MC simulations following the Metropolis scheme were used . 
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The total energy of interaction of the system can be written as 



tot 



Uhs + Ucoul + UpENE + UlJ + Uydw-i 



(1) 



i,i<j 



Uhsir) 



(2) 



where all the contributions of the pair potentials in Eq. (HI) are described in detail below. 

n 

Excluded volume interactions are modeled via a hard sphere potential Uhs defined as 
follows 

oo, for r < rcut 
0, for r > Tcut 

where r^ut = for the microion-microion excluded volume interaction, and r^-ut = a + a/2 
for the macroion-microion excluded volume interaction. Hence the center-center distance of 
closest approach between the macroion and a microion isro = a + o"/2 = 5o". 

The pair electrostatic interaction between two ions i and j (where i and j can be either 
a microion or the macroion) reads 



Ucoui{rij) = ±kBTlB 



ZiZj 



(3) 



where +(-) applies to charges of the same (opposite) sign and Ib = e"^ /AneoerkBT is the 
Bjerrum length corresponding to the distance at which two elementary charges interact 
with fc^T. To link our simulation parameters to experimental units and room temperature 
(T = 298K) we choose a = 4.25 A leading to the Bjerrum length of water Ib = 1.68a = 7.14 
A and to a macroion surface charge density of 0.14 ■ Cm~^. 

The polyelectrolyte chain connectivity is modeled by using a standard FENE potential 
in good solvent (see, e.g., which reads 



-nRn In 



UpENEir) = < 



1-- 

R? 



J 



oo, 



for r < Rq 



for r > Rn 



(4) 



where we chose k, = 27kBT/a'^ and Rq = 1.5a. The excluded volume interaction between 
chain monomers is taken into account via a purely repulsive Lennard- Jones (LJ) potential 
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given by 



UlAt) 



4e 



a 
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+ e, for r < 2^^ 



(5) 



0, for r > 2V6(j 

where e = /cbT. These parameters lead to an equihbrium bond length I — 0.98(7. 

An important interaction in the multilayering process addressed in this study is the non- 
electrostatic short ranged attraction, Uvdw, between the macroion and the PC chain. To 
account for this kind of interaction, we choose without loss of generality a van dcr Waals 
(VDW) potential of interaction between the macroion and a PC monomer that is given by 

6 



Uvdw{r) -ex 



vdw 



for r > ro, 



(6) 



r - ro + (7 

where Xvdw is a positive dimensionless parameter describing the strength of the attraction. 
Thereby, at contact (i.e., r = ro), the magnitude of the attraction is Xvdw^ = XvdwksT, and 
for Xvdw — 1, one recovers the standard attractive component of the LJ-potential [see Eq. 



TABLE I: Model simulation parameters with some fixed values. 



Parameters 



T = 298K 

Zm = ^0 

Z = 1 

a = 4.25 A 

Ib = 1.680- = 7.14 A 

a = 4.5<T 

ro = o + I = 5(7 

R = 60a 

N+ 

A_ 

NpE = N+ + N_ 

Nm 
Xvdw 



room temperature 
macroion valence 
microion valence 
microion diameter 
Bjerrum length 
macroion radius 
macroion-microion distance of closest approach 
radius of the outer simulation cell 
number of PCs 
number of PAs 
total number of polyelectrolyte chains 
number of monomers per chain 
strength of the specific van der Waals attraction 
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(0)]. Since it is not straightforward to link this strength of adsorption directly to experi- 
mental values, we therefore investigated different possible strengths of attraction, which are 
known to be realistic for soft matter systems. 

All the simulation parameters are gathered in Table HI The set of simulated systems 
can be found in Table IHI Single-particle moves were considered with an acceptance ratio 
of 20 — 30% for the monomers and 50% for the counterions. At equilibrium, the (average) 
length of the trial moves Ar are about 30cr for the counterions and O.lcr for the monomers. 
About 10^ to 10^ MC steps per particle were required for equilibration, and about 2 x 10^ 
subsequent MC steps were used to perform measurements. 

III. TARGET QUANTITIES 

Before presenting the results, we briefly describe the different observables that are going 
to be measured. Of first importance, we compute the radial density of monomers n±{r) 
around the spherical macroion normalized as follows 

R 

4:Trr^n±{r)dr = N±Nm (7) 



where (+)— applies to PCs (PAs). This quantity is of special interest to characterize the 
degree of ordering in the vicinity of the macroion surface. 



TABLE II: System parameters. The number of counterions (cations and anions) ensuring the 
overall electroneutrality of the system is not indicated. 



System 


NpE 


N+ 


A_ 




A 


1 


1 





80 


B 


2 


1 


1 


80 


C 


3 


2 


1 


80 


D 


4 


2 


2 


80 


E 


5 


3 


2 


80 


F 


6 


3 


3 


80 


G 


12 


6 


6 


80 


H 


40 


20 


20 


10 
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The total number of accumulated monomers A^±(r) within a distance r from the macroion 
center is then given by 



where (+)— applies to PCs (PAs). This observable will be used for the study of the adsorp- 
tion of (i) a single PC chain (Sec. IIV|) and (ii) two oppositely charged polyelectrolytes (Sec. 



Another quantity of special interest is the global net fluid charge Q{r) which is defined 
as follows 



where n_|_ (n_) include the density of all the positive (negative) microions. Thus, Q{r) 
corresponds to the total fluid charge (omitting the macroion bare charge Zm) within a 
distance r from the macroion center, and at the cell wall Q{r = R) = Zm- Up to a factor 
proportional to [Q{r) — Zm] gives (by simple application of the Gauss theorem) the 

mean electric field at r. Therefore Q{r) can measure the strength of the macroion charge 
screening by the charged species present in its surrounding solution. 

IV. SINGLE CHAIN ADSORPTION 

In this part, we study the adsorption of a single long PC chain (system A) for different 
couplings Xvdw Experimentally this would correspond to the process of the first polyelec- 
trolyte layer formation. 

The monomer density n+(r) and fraction Nj^{r)/Nm are depicted in Fig. H^a) and Fig. 
Q]^b), respectively. The density ?T,+ (r) near contact (r ~ ro) increases considerably with Xvdw 
as expected. At a radial distance of l.Sa from the macroion surface (i.e., r = ro + a), more 
than 97% of the monomers are adsorbed for sufficiently large Xvdw (here Xvdw > 3) against 
only 78% for Xvdw = 0. 

The net fluid charge Q{r) is reported in Fig. |21 In all cases we observe a macroion charge 
reversal (i.e., Q{r)/ZM > 1), as expected from previous studies j^, addressing only 
Xvdw = 0. The position r = r* at which Q{r*) gets its maximal value decreases with Xvdw, 
due to the Xvdw-Gnh.a.iaced adsorption of the chain. This overcharging increases with Xvdw, 
since the gain in energy by macroion-monomer VDW interactions can better overcome (the 
higher Xvdw) the cost of the self-energy stemming from the adsorbed excess charge. Note 
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FIG. 1: Monomer adsorption profiles of a single PC-chain (system ^4) at different Xvdw couplings, 
(a) radial density n+(r). (b) fraction of monomers N^(r)/Nm- 

that the maximal value of charge reversal of 100% allowed by the total PC charge (i.e., 
Q{r*)/ZM = 2) can not be reached due to a slight accumulation of microanions. 

Typical equilibrium configurations can be found in Fig. El For all Xvdw values, there 
is a wrapping of the chain around the macroion. In parallel, one can clearly see that the 
formation of chain loops is gradually inhibited by increasing Xvdw 

Although all the obtained results are intuitively easy to understand, they will turn out 
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(r-r^)/G 

FIG. 2: Net fluid charge Q{r) for system A at different Xvdw couplings. The horizontal line 
corresponds to the isoelectric point. 
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FIG. 3: Typical equilibrium configurations for a single PC adsorbed onto an oppositely charged 
macroion (system A). The little ions are omitted for clarity, (a) Xvdw = (b) Xvdw = 1 (c) 
Xvdw = 2 (d) Xvdw = 3 (e) Xvdw = 5. 

to be helpful in order to have a quantitative analysis of the effect of an extra short-range 
attraction already on the level of a single chain adsorption. 

V. ADSORPTION OF TWO OPPOSITELY CHARGED POLYELECTROLYTES 

We now consider the case where we have additionally a PA chain (system B), so that we 
have a neutral polyelectrolyte complex (i.e., one PC and one PA). Experimentally this would 
correspond to the process of the second polyelectrol5d;e layer formation (with system A as 
the initial state). We stress the fact that this process is fully reversible for the parameters 
investigated in our present study. In particular, we checked that the same final equilibrium 
configuration is obtained either by (i) starting from system A and then adding a PA or 
(ii) starting with no chains and then adding the two oppositely charged polyelectrolytes, 
together with their counterions, simultaneously. 
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FIG. 4: Monomer adsorption profiles of two polyelectrolyte chains (system B) at different Xvdw 
couplings. The solid and dashed lines correspond to PC and PA monomers, respectively, (a) radial 
density n±{r). (b) fraction of monomers N±{r)/Nm- 

The monomer density ?T,-|-(r) and fraction N±{r)/Nm are depicted in Fig. Ufa) and Fig. 
Efb), respectively. The corresponding microstructures are sketched in Fig. El The density 
of PC monomers n+(r) near contact (r ~ tq) increases considerably with Xvdw as expected. 
However, a comparison with system A (see Fig. indicates that the adsorption of PC 
monomers (at given Xvdw) is weaker when an additional PA is present. This is consistent 
with the idea that the PC chain tends to build up a globular state by getting complexed 
to the PA chain. This feature is well illustrated in Fig. El More precisely, for sufficiently 
small Xvdw ^ 1, the polyelectrolyte globular state is highly favorable compared to the "flat" 
bilayer state (see also Fig. Nevertheless, at sufficiently large Xvdw ^ 2, the first layer 
made up of PC monomers is sufficiently stable to produce a second layer made up of PA 
monomers. Thereby, the two chains wrap around the macroion. As far as the PA monomer 
adsorption is concerned. Fig. jH shows that n_(r) always increases with Xvdw- For Xvdw = 0, 
the polyelectrolyte complex is very close to the globular polyelectrolyte hulk state (i.e., in 
the absence of the macroion). This a non-trivial result, since naively one would expect a 
"true" multilayering for any Xvdw- It is only for large Xvdw ^ 3 that one gets a true bilayer 
formation, where there is a pronounced peak in n^{r) around r — tq ~ o". 

It is useful to introduce the following dimensionless interaction parameters Tm = ZmZ^-^, 
which measures the strength of the macroion-PC electrostatic attraction, and = 2'^^ 
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FIG. 5: Typical equilibrium configurations for one PC (in white) and one PA (in red) adsorbed 
onto the charged macroion (system B). The little ions are omitted for clarity, (a) Xvdw = (b) 
Xvdw = 1 (c) Xvdw = 2 (d) Xvdw = 3 (e) Xvdw = 5. 

which controls the PC-PA complex interaction. For large values of F^ the bulk complex will 
always be in a globular state, since then the Coulomb interaction will give rise to a chain 
collapse, similar to those seen in polyelectrolyte systems. Thus, for a sufficiently large value 
of Fto/Fm at given Xvdw, the globular state will always be favorable and no bilayering can 
occurs. In this case unwrapping occurs, similarly to the microstructures depicted in Fig. 
Efa) and Fig. El^b). On the other hand, we find at fixed parameters F^ and Tm, that one 
needs a s sufficiently large value xldw^ order to achieve bilayering. 
One can summarize these important results as follows: 

• The equilibrium bilayering process on a spherical charged colloid with long polyelec- 
trolyte chains requires a sufficiently strong extra short-ranged macroion-PC attraction. 

A closer look on Fig. Efb) reveals a further non-trivial behavior in the profiles of N±{r) 
at high Xvdw- Very close to the macroion surface we always have a monotonic behavior of 
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the fraction of adsorbed PC [N^{r;Xvdw)] and PA monomers [N^{r]Xvdw)] with respect to 
Xvdw as it should be. However, for a certain distance away from the surface we observe an 
unexpected behavior where iV+(r; Xvdw = 3) > iV+(r; Xvdw = 5) as well as N_{r; Xvdw = 3) > 
N_{r; Xvdw = 5). One can qualitatively explain this effect by the onset of the formation of 
one (or several) polyelectrolyte microglobule(s), i.e., small cluster (s) of oppositely charged 
monomers [see Figs. |Sl(d) and (e)]. This is indeed possible because at high Xvdw in principle 
more PC ( and consequently also PA) monomers want to get close to the macroion surface. 
Already for neutral chains a two dimensional flat adsorbed chain needs a high (surface) 
binding energy. Compared to bulk conformations the chain entropy is roughly reduced by 
kBTN\\i{qct=2/ (ld=?^- Here q is the effective number of conformational degrees of freedom 
of a bond, giving that \ii.{qd=2 / (ld=3) = 0{1). Thus local microglobules that induce a small 
local desorption, are entropically much more favorable. Also, on the level of the energy, an 
increase of q concomitantly favors the PC-PA microglobule. 

The net fluid charge Q{r) is reported in FigEl For Xvdw ^ 2 the macroion gets even 
overcharged and undercharged as one gets away from its surface, whereas for low Xvdw no 
local overcharging occurs. Again, at high Xvdw the strength of the charge oscillation is not a 
monotonic function of Xvdw where we observe a higher local overcharging (and undercharging) 
with Xvdw = 3 than with Xvdw = 5. This latter feature is fully consistent with the profiles 
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(r-r^)/o 

FIG. 6: Net fluid charge Q{r) for system B at different Xvdw couplings. The horizontal line 
corresponds to the isoelectric point. 
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of N±{r) [see Fig. Efb)] previously discussed However this onset of local (surface) 

microglobules (for Xvdw = 5) is not strong enough to produce a non-monotonic behavior of 
r* with respect to Xvdw 



VI. MULTILAYERING 



We now turn to the case where there are many polyelectrolytes (with NpE > 3) in the 
system. We recall that when Xvdw 7^ 0, the VDW interaction concern monomers of all PCs 
lying within the range of interaction. To keep the number of plots manageable, we will 
present our results obtained for Xvdw = and Xvdw = 3. The case Xvdw = is (conceptually) 
important since it corresponds to the situation where only purely electrostatic interactions 
are present. The other case Xvdw = 3 seems to be a reasonable choice, since we found 
a "true" bilayering for that value. Moreover, such a strength should be easily accessible 
experiment ally. 



x = o 




FIG. 7: Radial niononier density for the systems .A — G with \vdw 

= 0. The solid and dashed 

lines correspond to n+(r) and n-{r), respectively. The number of PC and PA chains is indicated. 
The plots of n±{r) for the systems A (IPC) and B (IPC-lPA) are again reported here for direct 
comparison. 
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FIG. 8: Typical equilibrium configurations for many polyelectrolyte chains adsorbed onto the 
charged macroion with Xvdw = 0. The PC monomers are in white and PA ones in red. The 
little ions are omitted for clarity, (a) 2PC-1PA (system C) (b) 2PC-2PA (system D) (c) 3PC-2PA 
(system E) (d) 3PC-3PA (system F) (e) 6PC-6PA (system G). 

A. Adsorption with Xvdw = 

The density profiles of n±{r) for the systems A — G (with Xvdw = 0) are reported in Fig. 
[7|and the corresponding microstructures are sketched in Fig. |H| 
Figure [7| shows that when the total polyelectrolyte charge, 

QPE = - N_)M^e, (10) 

is zero, the density of PC monomers n+(r) near contact is lower than when charge Qpe = 
N^e (recalling that our systems are such that Qpe = or N^e). This tendency [lower n_|_(r) 
near contact with Qpe = 0] gradually decreases as the total number NpE of polyelectrolytes 
is increased. In particular for the system G where NpE = 12 and Qpe = 0, the density 
?T,_|_(r) near contact is nearly identical to that of systems A, G and E where Qpe = N^g. 
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On the other hand, when Qpe = then n+(r) near contact is nearly independent of 

NpE- The height of the peak in the PA monomer density n_(r) increases monotonically 
with NpE- Concomitantly, a third layer made of PC monomers builds up for high enough 
NpE- This multilayering is especially remarkable for NpE = 12 (system G). 

All these features concerning the first layer structure can be rationalized with simple 
ideas. When Qpe = 0, then the resulting global attraction between the macroion and the 
polyelectrolyte complex is much weaker than when Qpe = NmG. In this latter situation 
where Qpe = Nj^e, this excess charge carried by a PC chain leads to a relatively strong PC 
adsorption near the surface. The equilibrium configurations sketched in Fig. |H1 suggest a 
wrapping from the PC monomers when Qpe = N^e [see Fig. IHI^a) and (c)] and a (partial) 
unwrapping when Qpe = [see Fig. IHb), (d)]. Even for high NpE = 12 [see Fig. Et^e)], 
we can see this tendency of unwrapping leading to a polyelectrolyte-complex globular state. 
However, for symmetry reasons, when the total number of monomers is large enough as it is 
the case with NpE = 12, the distribution of the polyelectrolyte complex gets more isotropic 
leading to a weaker unwrapping at Qpe = 0. 

The collapse into a globular polyelectrolyte complex becomes even more spectacular when 
a is reduced (i.e., larger F^) |30|. In that case (not reported here), we found a wrapping 
(for Qpe = Nmc) similar to that depicted in Fig. |Hl^a) and (c), and a strong unwrapping 
(for Qpe = 0) where the compact neutral polyelectrolyte complex is adsorbed onto a small 
area of the colloid. 

The net fluid charge Q{r) is reported in Fig. IHl As expected one detects an overcharging 
and undercharging for Qpe = NmC and Qpe = 0, respectively. For Qpe = 0, the macroion 
is also locally overcharged very close to the macroion surface and its strength increases with 
NpE- On the other hand, the strength of the undercharging (occurring at the largest radial 
position r* of the extrema) at Qpe = is nearly independent of NpE- In parallel, the 
strength of the overcharging (occurring at the largest radial position r* of the extrema) 
measured at Qpe = NmC is also nearly independent of NpE (systems C and E). Moreover, 
our simulations show that for NpE > 2 the strength of the overcharging (with Qpe = Nmc) 
and undercharging (with Qpe = 0) have nearly the same amplitude, in qualitative agreement 
with experimental data. 
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B. Adsorption with Xvdw 7^ 

In this part, we consider the additional attractive VDW macroion-PC monomer interac- 
tion with Xvdw = 3. The same investigation as with Xvdw = is carried out here. 

The density profiles of n±{r) for the systems A — G (with Xvdw = 3) are reported in Fig. 
^Jand the corresponding microstructures are sketched in Fig. ITTl 

Figure ITUl shows that the density n^{r) near contact (for a given system) is about six 
times larger than that obtained at Xvdw = (compare Fig. [7j). When Qpe = NmC, the 
density n^{r) at contact (slightly) increases monotonically with NpE in contrast to what 
happened at Xvdw = where it was nearly independent of NpE- When Qpe = 0, we remark 
that the density n^{r) near contact is nearly independent of NpE (for NpE > 2) in contrast 
to what happened at Xvdw = 0. 

As far as the PA density n_{r) is concerned, the height of the first peak (for a given 
system) is about twice larger than that obtained at Xvdw = 0. This height is a monotonic 
function of NpE within a given regime of Qpe (here, either or Nme). Nevertheless, in 
general this height exhibits a non-trivial dependence on NpE, in contrast to our results with 
Xvdw = 0. For the systems B and C both containing a single PA chain (A^_ = 1), the height 
of the first peak in n_(r) is smaller with iV_|_ = 2 (system C) than with iV_|_ = 1 (system 

z = o 

2 



1.5 
0.5 



5 10 15 

(r-r^)/o 

FIG. 9: Net fluid charge Q{r) for the systems A — G with Xvdw = 0. The number of PC and PA 
chains is indicated. The plots for the systems A (IPC) and B (IPC-lPA) are again reported here 
for direct comparison. The horizontal hne corresponds to the isoelectric point. 
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B). This is again due to the formation of clusters of oppositely charged monomers that 
takes place above the first layer. This effect is more pronounced when the amount of PC 
monomers (at given N_) is larger (system C), leading to a local desorption of PA monomer. 
Those features are remarkable by a comparison of the snapshots of the systems B and C 
depicted in Fig. ]^d) and Fig. ITlT a). respectively. Similar arguments can be used for the 
systems D and E, where the same effect is found. At large NpE, the height of the first peak 
in n_(r) saturates as expected. 

For 3 < NpE < 6, our simulation shows that the formation of the third layer [i.e., the 
second peak in n-^{r) at r — tq ~ 2.6cr] is enhanced when Qpe = at fixed iV+. This 
effect can again be explained in terms of polyelectrolyte (micro)globules formation. Indeed, 
above the second layer, the formation of clusters made up of oppositely charged monomers 
is enhanced when the polyelectrolyte complex (seen by the underneath bilayer) is uncharged 
which corresponds to a state of charge Qpe = 0. 

It is interesting to see that with NpE = 12 one even gets a second peak (and not a 
shoulder) in n-{r), which is the signature of a fourth layer. This qualitatively contrasts 
with our findings at Xvdw = 0. Therefore, we conclude that the effect of an extra short- 
ranged macroion-PC attraction is crucial for the multilayering process. 

On a more qualitative level, it is very insightful to compare the microstructures obtained 
with purely electrostatic interactions {Xvdw = 0) sketched in Fig. |H1 with those obtained 
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FIG. 10: Same as Fig. [7|but with Xvdw = 3. 
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FIG. 11: Same as Fig. |Hlbut with Xvdw = 3. (a) 2PC-1PA (system C) (b) 2PC-2PA (system D) 
(c) 3PC-2PA (system E) (d) 3PC-3PA (system F) (e) 6PC-6PA (system G). 

with a short-ranged VDW macroion-PC interaction {xvdw = 3) sketched in Fig. El From 
such a visual inspection, it is clear that in all cases the adsorbed polyelectrolyte complex 
is flatter at Xvdw = 3 than at Xvdw = 0. An other important qualitative difference, is that 
the unwrapping occurring at Xvdw = with Qpe = [see Fig. |Hl^b) and (d)] is no longer 
effective when Xvdw = 3 [see Fig. ITlT b) and (d)]. In the same spirit, for a large number 
of chains {NpE = 12), the macroion surface is only partially covered by the PC monomers 
where some (large) holes appear [see Fig. IHe)], in contrast to what occurs at Xvdw = 3, 
where all the macroion surface is covered [see Fig. ITTT e)]. 

The net fluid charge Q{r) is reported in Fig. ^1 As expected one finds an overcharging 
and undercharging for Qpe = and Qpe = 0, respectively. Now one can get a local 
overcharging larger than 100% (i.e., Q{r)/ZM > 2) due to the VDW attraction that can lead 
to a first layer with many PC chains. For systems C and E we see that the overcharging at 
the third layer is around 50% and nearly independent of NpE- 

On the other hand, the strength of the undercharging (occurring at the largest radial 
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FIG. 12: Same as Fig. 1^1 but with Xvdw = 3. 

position r* of the extrema) at Qp^; = decreases with increasing NpE, providing a gradually 
weaker driving force for the subsequent adsorption of the PC chain. On the basis of our 
results with Xvdw = 0, we expect that the strength of the undercharging at Xvdw = 3 (for 
larger Npe) should stabilize around 50%. So it appears that the oscillation of under- and 
overcharging are not 100%, but rather close to 50%. This is probably sensitive to the specific 
model parameters chosen. What can be stated from our data is, that there is no reason to 
find a generally applicable overcharging fraction. Especially for the case of relatively small 
colloids, the results will strongly depend on the specific system parameters, which are both 
of electrostatic as well as non-electrostatic nature. 

VII. CASE OF SHORT CHAINS 

We now investigate the effect of chain length dependence. In this case the adsorption of a 
single chain does not necessarily produce an overcharging since the chain length {Nm = 10 - 
system H) is too short. The density profiles of n±{r) are reported in Fig. HSlfor various Xvdw, 
and the corresponding microstructures are sketched in Fig. EI In the purely electrostatic 
regime {Xvdw = 0), the polyelectrolyte adsorption is weak and it significantly increases with 
Xvdw However, for all reported cases, we only observe a bilayering in contrast with previous 
long chain systems (compare Fig. ^1 with Fig. [7| and Fig. ITUI) where thereby a true 
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Short chains 




FIG. 13: Radial nionomer density for short polyelectrolyte chains (system i?) at different x.vdw 
couphngs. The black and gray lines correspond to PC and PA monomers, respectively. 

multilayering was reported. 

In addition we observe several globally neutral polyelectrolyte complexes in the bulk, 
whose number decreases with Xvdw (see Fig. IT^ . This feature was inhibited for long chains 
due to the strong PC-PA binding energy that keeps all the chains near the macroion surface. 
At sufficiently strong Xvdw [see Fig. IT^ c) with Xvdw = 5], the macroion area gets largely 
(and uniformly) covered by the PC chains, leading to a strong bilayering. Nevertheless, due 




FIG. 14: Typical equilibrium configurations for short PC (in white) and PA (in red) chains 
adsorbed onto the charged macroion (system H). The little ions are omitted for clarity, (a) 

Xvdw = (b) Xvdw = 3 (c) Xvdw = 5 
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Short chains 




FIG. 15: Net fluid charge Q{r) for system H at different Xvdw couplings. The horizontal line 
corresponds to the isoelectric point. 

to the weak PC-PA binding energy, the formation of additional layer seems to be prohibited 
in contrasts to what was observed at Xvdw = 3 with systems D and E that contain a 
similar number of monomers. Those observations lead us to the relevant conclusion that 
multilayering with short chains (if experimentally observed on a charged colloidal sphere) 
must involve additional non-trivial driving forces like specific PC monomer-PA monomer 
interactions that are not captured by our model. This again seems to be in agreement with 
the arguments presented in Ref.j?! which argue against a stable thermodynamic equilibrium 
complex when there is excess polyelectrolytes present. 

The net fluid charge plotted in Fig. indicates that only one charge oscillation (around 
the isoelectric point) is obtained in contrast to what can happen with longer chains. Again, 
here the driving force for the bilayer formation is the overcharging that increases with Xvdw 



VIII. CONCLUDING REMARKS 



We have carried out MC simulations to study the basic mechanisms involved in forming 
equilibrium polyelectrolyte complexes on a charged colloidal sphere. This work emphasizes 
the role of the short-range Van der Waals-like attraction (characterized here by Xvdw) between 
the spherical macroion surface and the oppositely charged adsorbed chain(s). 

It was demonstrated that, for the bilayering process involving two long oppositely charged 
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chains at fixed and Tm, it is necessary to have a sufficiently high Xvdw In particular, 
below a certain value of Xvdw = xldw^ ^ dense adsorbed polyelectrolyte globule is obtained, 
whereas above xldw ^ bilayer builds up. 

The same qualitatively applies to the case of many (more than two) long polyelectrolytes. 
In a purely electrostatic regime (i.e., Xvdw = 0) one can never obtain a true (uniform) 
multilayering. However, by increasing Xvdw, one gradually increases the polyelectrolyte 
(polycation and polyanion) chain adsorption ultimately leading to a true multilayering where 
the macroion is uniformly covered. Nonetheless, at given Xvdw and especially for small Xvdw, 
the polyelectrolyte globular state is always favored when its net charge is zero. 

As far as the short chain case is concerned, it was shown that even bilayering can not be 
reached within the pure electrostatic regime. Only at higher Xvdw (higher than those coming 
into play with long chains), one recovers a bilayer formation. However, multilayering (beyond 
bilayering) with very short chains seems to be very unlikely within our model. The large 
complex would not be thermodynamically stable and dissolve into smaller charge neutral 
polyelectrolyte complexes, consistent with the ideas presented in Ref. [3]. 

As an overall conclusion, our results clearly demonstrated that besides an overcharging 
driving force [i.e., successive macroion (effective) charge reversal by successive polymer lay- 
ering), the stability of the polyelectrolyte multilayer is strongly inffuenced by the speciffc 
macroion-polyelectrolyte short range attraction. This statement should at least hold for the 
investigated cases of equilibrium structures. 

A future study should systematically study other important effects, such as chain ffexibil- 
ity, specific interchain monomer-monomer interaction, microions valency etc... Nevertheless, 
our present findings hopefully will generate some further systematic studies to explore the 
effects of non-electrostatic effects for the layer-by- layer deposition technique. 
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